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ABSTRACT 

We perform a detailed X-ray study of the filaments surrounding the brightest cluster galaxies 
in a sample of nearby galaxy clusters using deep Chandra observations, namely the Perseus, 
Centaurus and Virgo clusters, and Abell 1795. We compare the X-ray properties and spectra of 
the filaments in all of these systems, and find that their Chandra X-ray spectra are all broadly 
consistent with an absorbed two temperature thermal model, with temperature components 
at 0.75 and 1.7 keV. We find that it is also possible to model the Chandra ACIS filament 
spectra with a charge exchange model provided a thermal component is also present, and the 
abundance of oxygen is suppressed relative to the abundance of Fe. In this model, charge ex¬ 
change provides the dominant contribution to the spectrum in the 0.5-1.0 keV band. However, 
when we study the high spectral resolution RGS spectrum of the filamentary plume seen in 
X-rays in Centaurus, the opposite appears to be the case. The properties of the filaments in 
our sample of clusters are also compared to the X-ray tails of galaxies in the Coma cluster and 
Abell 3627. In the Perseus cluster, we search for signs of absorption by a prominent region of 
molecular gas in the filamentary structure around NGC 1275. We do find a decrement in the 
X-ray spectrum below 2 keV, indicative of absorption. However the spectral shape is inconsis¬ 
tent with this decrement being caused by simply adding an additional absorbing component. 
We find that the spectrum can be well fit (with physically sensible parameters) with a model 
that includes both absorption by molecular gas and X-ray emission from the filament, which 
partially counteracts the absorption. 
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1 INTRODUCTION 

Filaments of material which are bright in their X-ray and Ha emis¬ 
sion have been found surrounding the brightest cluster galaxies 
(BCGs) in a number of galaxy clusters. Detailed studies of indi¬ 
vidual, nearby systems such as NCG 1275 in the Perseus cluster 
have revealed that the histories of these filaments are intertwined 
with those of the X-ray cavities in the intracluster medium (ICM) 
produ ced from AGN feedback (e.g. lFabian et al.lE003l . lHatch et al.l 
2006). The filaments are dragged outwards from the centre of the 
BCG by the buoyantly rising bubbles, before falling back. The mor¬ 
phologies of these filaments can be used as streamlines, allowing us 
to trace the flow of gas in the ICM. 

The remarkable linearity of the filaments around NGC 1275, 
some of which are over 6kpc long yet only 70 pc wide, indicates 
that they are supported by magnetic fields which stabilise them 
again st tidal shear and dissipation into the surrounding, hotter ICM 
dFabian et al .I|200f3 ). 


* Email: swalker@ast.cam.ac.uk 


The excitation mechanism of the filaments remains uncertain. 
Around NGC 1275 in the Perseus cluster, photoionisation is an 
unlikely mechanism as the sources of ionizing radiation are not 
sufficiently luminous, and there are no obvious clusters of stars 
associated with the filaments. The filaments also have an emis¬ 
sion spectrum which is different from any known Galactic neb¬ 
ula. Around other BCGs, some filaments are seen to be associ¬ 
ated with ongoing star formation (for example, ; McNamara~etalJ 


2004 1 0’ Pea et al 1 2004 lo ’ Pea et all2010LlMcDonald et all201lL 

Tremblay et all2015l andlDonahue et al.ll2015h. 


Ijohnstone et all d2007h studied the H 2 emission from the fila¬ 
ments around NGC 1275 using Spitzer and ground-based observa¬ 
tions and found that the level excitation energy was correlated with 
the derived excitation temperature. It was then found that this cor¬ 
relation co uld be accounted for by collisional heat ing from ionizing 
particles dFerland et al.ll2008llFerland et al.ll2009it . with agreement 
within a factor of 2 over a wide range of line ratios. It is most likely 
that the optical/infrared/radio emission-line spectrum is produced 
by gas exposed to ionizing particles. 


iFabian et alj d201 lh studied a prominent filament to the north 
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Figure 1 . Images of the filaments (in H a except for A1795, for which the Chandra X-ray data are shown) and the regions used to study them. Green boxes 
show the ‘on’-filament regions, while the cyan boxes for the local background (‘off'-filament) regions used for the X-ray analysis of each filament. The clusters 
shown are Perseus (top left), Abell 1795 (top right), Centaurus (bottom left) and Virgo (bottom right). 


of NGC 1275 and found that the surface radiative flux from the fil¬ 
ament is comparable to the energy flux on itjrom the surrounding 
hot gas. Given this remarkable similarity, iFabian et alj j 201 ll) ar¬ 
gued that the most likely energy source of these i onizin g particles 
is from the hot gas surrounding the filaments. [Fabian et al.1 1201 ll) 
proposed that the hot gas surrounding the cold filaments is able 
to pe netrate into them through the process of reconnection diffu¬ 
sion lLazarian et alJl20l(il . lLazarian et al.ll201lh . which provides a 
mechanism through which the hot gas can penetrate into the cold 
filaments despite the magnetic fields which support the filaments. 
This penetration process provides the secondary electrons within 
the cold filaments which excite the observed submillimetre radia¬ 
tion though UV emission. The interpenetration of the cold gas in 


the filaments by the surrounding hot gas provides another mecha¬ 
nism through which the hot gas can lose energy, which may mean 
that gas can cool from the hot ICM at a rate higher than that de¬ 
duced from just X-ray spectroscopy. 

In this mechanism, the hot ionised gas surrounding the fila¬ 
ments penetrates and mixes with the cold neutral gas in the fila¬ 
ments. Charge exchange will then dominate, with electrons being 
transferred from the cold neutral gas into excited levels of the im¬ 
pinging hot ions. X-ray line emission then results from the cascade 
of the electrons down to lower energy levels. This process has been 
observed in other astrophysical environments, such as for exam¬ 
ple the X-ray emission seen from comets caused by their hydrogen 
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Table 1. Coordinates, redshift and absoiption data of objects used in analy¬ 
sis. 


Object 

RA 

Dec 

z 

N H /(cm- 2 ) 

Perseus 

03hl9m47.2s 

+41d30m47s 

0.017900 

1.4xl0 21 

Centaurus 

12h48m47.9s 

-41dl8m28s 

0.011400 

8.30xl0 20 

M87 

12h30m49.4s 

+12d23m28s 

0.004283 

2.07 xlO 20 

A1795 

13h48m53.0s 

+26d35m44s 

0.062476 

1.21 xlO 20 

Coma 

12h59m48.7s 

+27d58m50s 

0.023100 

8.43xl0 19 

A3627 

16hl4m22.5s 

-60d52m07s 

0.015700 

1.81X10 21 


haloes interacting with highly charged ions in the solar wind (e.g. 
iDennerl et alii 1997 )._ 

[Fabian et alj 1201 ll) found that X-ray spectrum of the north¬ 
ern filament of NGC 1275 was dominated by line emission, and 
could be well fitted with a two temperature model with compo¬ 
nents at ~0.7 and ~ 1.5 keV. In the proposed reconnection diffusion 
model, charge exchange is expected to contribute to the soft X-ray 
emission of the filaments through Fe-L and O emission lines. How¬ 
ever, [F^ian^r^] |2oTl!) found that the predicted soft X-ray flux is 
around a factor of 30 lower than that which is actually observed, if 
it is assumed that each incoming ion has only one chance of charge 
exchange per stage. One possibility is that the cold filaments are 
full of tiny strands, which could increase the rate of interaction re¬ 
sulting in charge exchange emission. 

Here we explore the X-ray spectra of the filaments around the 
BCGs of nearby galaxy clusters for which very deep, high spatial 
resolution Chandra observations are available. These are namely 
the Perseus cluster, the Centaurus cluster, the Virgo cluster and 
Abell 1795. We investigate whether or not the X-ray spectra of the 
filaments in these different systems are similar, to try understand 
and constrain the underlying mechanism responsible for their soft 
X-ray emission. We test whether the latest charge exchange model 
ACX (AtomDB Charge Exchange model. [Smith et al.l2012h is able 
to to model the X-ray spectra of the filaments in Perseus, and the 
RGS spectrum of the filamentary plume in the Centaurus cluster. 

We use a standard ACDM cosmology with Ho = 70 km s -1 
Mpc -1 , Qm = 0.3, Ha=0.7. All errors unless otherwise stated are 
at the 1 a level. 


2 DATA REDUCTION 

We stacked all of the available Chandra archive data on the Perseus 
cluster (890 ks), the Centaurus cluster (760 ks), the Virgo cluster 
(500 ks) and Abell 1795 (700 ks). The observations used are tab¬ 
ulated in table IA1I in appendix [A] and these were reduced using 
CIAO version 4.7. The events files were reprocessed using the task 
CHANDRA_REPRO. Light curves in the 0.5-7.0 keV band were 
examined, and the routine LC_SIGMA_CLIP was used to remove 
any periods of flaring where the count rate differed from the mean 
by more than 2 a. The resulting clean exposure times are tabulated 
in table lAll 


3 METHODOLOGY 

The filaments around the BCGs in the systems we study are most 
easily seen through their Ha emission. We therefore choose the re¬ 
gions to study based on these Ha images (except for A1795, for 



Figure 2. X-ray images of the galaxy tails studies, along with the spectral 
extraction regions for the tails (green) and the local background regions 
used (cyan). The top panel shows the galaxy tail in the Coma cluster, while 
the middle and bottom panels show galaxy tails in Abell 3627. 


which the X-ray image is used), defining the ‘on’ and ‘off’ fil¬ 
ament regions shown in the panels in Fig |T] These panels show 
Perseus (top left, Ha WIYN telescope image from lConselice et ail 
l200 li). Abell 1795 (top right, Chandra 0.5-1.5 keV image), Cen- 
tauru s (botto m left, ESO Multi-Mode Instrument Ha image from 
ICrawford et ai] [2005 ) and Virgo (bottom right, Ha HST WFPC2 
image), which were used for extracting X-ray spectra from the 
Chandra data. The ‘off’-filament regions, shown in cyan were used 
to provide a local background and were subtracted from the ‘on’- 
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Figure 3. Plot of temperatures of filaments in the Perseus cluster. Blue 
points mark the temperature of the hotter component, red points represent 
the lower temperature component. 


filament regions, scaling for the difference in extraction size area 
between the foreground and the background. This allows the X-ray 
spectra of the excess emission from the filaments to be explored. 
For the Ha image of Virgo, the galaxy emission from the BCG has 
been subtracted to emphasise the filaments. 

We also explore three X-ray galaxy tails seen in the Coma 
cluster and Abell 3627, to see how their spectra compare to the 
filament spectra. The regions used in the spectral analysis of these 
tails are shown in the images in Fig. [2] 

For each spectral extraction region, the X-ray spectra were ex¬ 
tracted from each individual Chandra observation, together with the 
appropriate ARF and RMF calibration files for each observation us¬ 
ing the tools MKWARF and MKACISRMF respectively. The spec¬ 
tra from all of the observations of a particular region were then 
combined and weighted appropriately using the ARFs to produce 
the stacked spectrum for each region, while the individual ARFs 
and RMFs for each region were also combined with the correct 
weighting. Due to differences in the effective areas of the ACIS-S 
and ACIS-I detectors, the data from these two detectors were not 
added together. Instead, when both ACIS-S and ACIS-I data were 
used for the same target, we found a total ACIS-S spectrum and a 
total ACIS-I spectrum, and these two were fitted simultaneously in 
XSPEC. 

We fitted each local background subtracted filament 
region with an absorbed two component thermal model, 
PHABS(APEC-t-APEC), as has been employed in iFabian et al .1 
J20 1 ill . The metal abundances of the two components were tied 
together, allowing the temperatures and normalisations of each 
APEC component to be free. The column density was fixed to the 
LAB survey value dKalberla et al JI 2 OO 5 I) for each cluster studied, 
and these are tabulated in table Q] along with the values to which 
the redshift was fixed when fitting. 


3.1 Perseus 

For the Perseus cluster, the Chandra data are sufficiently deep and 
the filament emission is sufficiently bright that we can study each 



Energy (keV) 

Figure 4. The spectrum of all filaments (stacked) in Perseus that have been 
fitted with a two-temperature apec model. 


filament individually. As shown in the top left hand panel of Fig. [7] 
we identify 12 filaments in Perseus which we study individually. 

The best fitting parameters of the two thermal component 
model for the 12 filaments in Perseus are tabulated in table [2] and 
the best fitting temperatures are plotted in Fig. [3] In all of the sys¬ 
tems the two temperature model provides a good fit. The lower tem¬ 
perature components of the 12 filaments are all in good agreement, 
and are all reasonably consistent with the average temperature of 
0.77 keV (shown by the red line in Fig. [3}. The spread of the higher 
temperature component is higher, but again for all of the filaments 
these are in reasonable agreement with the average upper tempera¬ 
ture of 1.78 keV (shown by the blue line in Fig. [3}. 

Since the individual filaments appear to have similar spectral 
shapes, we stack them all together to produce a spectrum with the 
highest possible signal to noise, and this is shown in Fig.[4] showing 
components of the best fitting two component thermal model. 


3.2 Virgo, Centaurus and Abell 1795 filaments 

For the other clusters we study, the Chandra data are not sufficiently 
deep to allow individual filaments to be studied. Instead we obtain a 
stacked X-ray spectrum of all of the filaments in each cluster, and fit 
these with the same absorbed two thermal component model as be¬ 
fore. In each case the model provides a good fit to the data, and the 
best fitting temperature components are plotted in Fig. [5] In each 
case the lower temperature components appear reasonable consis¬ 
tent with one another. There is more scatter in the higher tempera¬ 
ture component, but again these are in broad agreement. 


3.3 Galaxy tails in the Coma cluster and Abell 3627 

The spectra of galaxy tails in the Coma cluster (originally pre¬ 


sented in 


Sanders et al. 2014h and Abell 3627 (originally 


pre¬ 


sented in ISun et al J l2Qlol and IZhang et al.1 l2013h were extracted 
as well, for comparison with the filaments in other clusters. The 
tail in the Coma cluster (Fig. [2] top panel), can be fitted with a 
two-temperature model, with a colder component of temperature 
0.221 qq 3 keV and hotter component of 1.22+o'lg keV. The abun¬ 
dance is 0.7 !q'q 5 Zq . For the first tail in Abell 3627 (Fig.[2] middle 
panel), we find that the best fit two temperature model is consistent 
with the findings of lSun et~al] l l2010l ). with components at 0.7(tg'J 
keV and lT^g’j keV and a best fit metal abundance of 0.2lg jZg. 
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Table 2. Fitting parameters for all filaments in the Perseus cluster. 


Number 

Temperature 1 (keV) 

Temperature 2 (keV) 

Abundance (solar) 

Norm 1 



Norm 2 



Norm ratio 

i 

0-83++ 

2 ifi+ 0 - 35 
z,1D -0.25 

0.43lg;£ 

^• DD —1.56 

X 

10" 6 

2 02"*”°' 34 
0.32 

X 

icr 5 

4 33+ 2 - 60 
^•oo_i 60 

2 

0-82±g;“ 

1 c:q+ 0-19 

-L.oO-o 24 

i.54lg;?8 

4 13+ 3 - 62 

H:.eo_2 69 

X 

to- 6 

8 01 +4 ‘ 91 
°* U1 _4.42 

X 

10~ 6 

1 94+ 2 - 08 
1,y ^-1.65 

3 

0-69±°;“ 

1 4Q+ 0 - 21 

1 -^ y -0.14 

1 f)0”*" 2 ‘ 44 

l.UU-o.54 

6.i9ifii 

X 

10" 6 

1 QQ+ 1 - 23 

i.yy_ 0- 7o 

X 

icr 5 

0 91 +3.66 
t3Zi —1.35 

4 

1-03++ 

2 3fi"*~ 1,22 
z,ou —0.59 

n 1O+0.15 

U-10 —0.07 

1.201J;8J 

X 

10- 5 

1 iq+0.46 
l-iy_0.64 

X 

icr 5 

0.99+® 2 

5 

o-8o ±g;ii 

1 40”*”°' 50 

0.23 

9 oc+2.30 
z * OJ —1.68 

1 13+ 119 
x -- LO —0.84 

X 

lO" 6 

2 04”*" 4 ' 27 
0.80 

X 

10" 6 

-f.8lll.51 

6 

n Qi TO. 1.4 
uyi —0.13 

1 QQ+ 5 - 51 
i.yy_074 

0-16++ 

1 10+117 
1 - 1 °-0.69 

X 

10" 5 

1 DrFO .62 

l.U8_ 0 95 

X 

io~ 5 

0.96+;“ 

7 

0.661°;“ 

l.DO-o 06 

0.69+“ 

1 07 + °' 21 
i - u/ -0.17 

X 

10- 5 

1 Dl” 1 ” 0 ' 14 

1-U1-0.13 

X 

10- 4 

9.44+; 20 3 

8 

0.66l°;i° 

1 65+ 0 - 08 

1 * oo -0.11 

0 49+ 016 

U -^ Z -0.13 

c 9q+1.31 
O.ZO-i 02 

X 

10" 6 

5 ni+°- 82 
°- ux —0.68 

X 

io~ 5 

9.581 2 ;® 7 

9 

rj 77+0.66 

U - ‘ ‘-0.09 

1 7Q+0.20 

' u —0.15 

0 44+ 0-20 
^•”—0.14 

1 1 9+0.35 
1,lz -0.29 

X 

10- 5 

q 71 +0.66 

' x -0.54 

X 

icr 5 

3 31+ 117 
o,o - l -0.97 

10 

rj 77+O.O6 

U -' 1 -0.06 

o 90+O.27 

0.20 

1 1 7+0.56 

1,1 ' —0.36 

6 ^7+ 2 - 28 

O.D/-1 78 

X 

lO" 6 

6 32 +1 ' 17 

o.oz_i 15 

X 

10- 5 

9 62"*” 3 ' 78 
y ‘ oz -3.14 

11 

0.86+ 22 

1 60” 1 ” 0,53 
i.OU-o 62 

2.18+;“ 

9 75 +7 ' 15 

1 °- 4.87 

X 

10- 7 

1 51+ 4 - 39 
lo± -1.03 

X 

10" 6 

1 fi c+4.66 
i,00 -1.34 

12 

o.69l°;“ 

1 09+O.24 
± - oz —0.15 

5.00++ 

1 14 - *" 1 ’ 22 
■ L,lf4 -0.22 

X 

lO" 6 

77+4.64 

1 1 -0.52 

X 

10" 6 

cc+4.40 

1.OO-0 54 
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Figure 5. Plot of temperatures of the best fit two temperature models of 
filaments in Perseus, M87, Centaurus and Abell 1795, and of the galaxy 
tails in Coma and Abell 3627. Blue points mark the temperature of the hotter 
component, red points represent the lower temperature component. 


The two temperature fit to the second tail in Abell 3627 is con¬ 
sistent with the findings of IZhang et alj d2013h , with temperature 
components at 0.5lo i6 keV and I. 2 I 015 keV and an abundance 
of0.3lg;? Z 0 . These temperature are also plotted and compared in 

Rg. a 


4 COMPARISON TO CHARGE EXCHANGE MODELS 

4.1 Comparing Chandra ACIS spectra of Perseus filaments 

Here we test whether the X-ray spectra of the filaments can be mod¬ 
elled with the charge exchange model ACjf] dSmith et ail 1 2012l ). 
The ACX model in the energy range of interest is shown in the left 
hand panel of Fig. [6] with the main emission lines labelled. The 
result of folding this model through the Chandra response is shown 

1 http://www.atomdb.org/CX/ 


in the right hand panel of Fig. [6] In this charge exchange model we 
have set the equilibrium ion population temperature to the temper¬ 
ature of the hot ICM surrounding the filaments, which is 4 keV. 

In Fig. [7] we compare the stacked spectrum of all of the fil¬ 
aments in Perseus to the shape of the charge exchange model 
ACX. To begin with we used the LAB survey column density of 
0.14xl0 22 cm -2 , for which the spectrum is shown as the solid 
black line in the top left hand panel of Fig. [7] We see that the shape 
of the charge exchange spectrum below 1.0 keV disagrees signifi¬ 
cantly with the observed X-ray spectrum of the stacked filaments, 
due to the OVII and OVIII lines being far brighter in the charge 
exchange model. 

We then allow the column density to be a free parameter, to see 
if this can suppress the OVII and OVIII lines sufficiently to produce 
an adequate fit to the observed spectrum. The resulting attempted 
fit is shown as the solid red line in the top left hand panel of Fig. [7] 
for which the column density is 0.36xlO 22 cm -2 . We see that the 
model remains a very poor fit to the data. 

Next we add an absorbed APEC component to the model, so 
that the total model is PHABS(APEC-l-ACX), to see if we can ob¬ 
tain an acceptable fit. The resulting fit is shown in the right hand 
panel of Fig. [7] where we have again allowed the absorption of 
the charge exchange model (dashed blue curve) to be a free pa¬ 
rameter. The fit improves on the addition of an APEC component 
but remains very poor, as can be seen from the residuals, with a re¬ 
duced X' 2 of 2.1 for 166 degrees of freedom in the 0.5-3.0 keV band. 
In this fit the column density is 0.2+q} x 10 22 cm -2 , while the 
APEC component, (which is the dominant component) has kT=1.4 
keV and an abundance of 0.2 Zq. 

Next we allow the relative abundances of O, Ne, Mg, Si 
and Fe to be free parameters, and try to fit the filament spectrum 
using an absorbed variable abundance charge exchange model, 
PHABS(VACX). The best fit is shown in the bottom left panel 
of Fig. [7] The best fit relative abundances for 0:Ne:Mg:Si:Fe are 
0.25:0.64:0.2:0.4:1.5. The enhancement of Fe relative to O and Ne 
dramatically improves the fit below 1 keV, and can correctly model 
the ‘step’ shaped feature between 0.7 and 1.0 keV. The best fit col¬ 
umn density is 0.231qqi x 10 22 cm -2 , and the quality of the fit 
is only slightly reduced if this is fixed to the LAB survey value of 
O. I47 q'qi x 10 22 cm -2 . However the fit above 1 keV remains poor, 
most notably the deficiency of the model in the 1-1.3 keV range, 
leading to a very poor goodness of fit statistic (reduced \ 2 =2.98 
for 102 degrees of freedom, fitting to the 0.5-2.0 keV band). 


© 0000 RAS. MNRAS 000, 000-000 

















6 S. A. Walker et al. 




Energy (keV) 

Figure 6. Left :The ACX spectral model, with the main lines in the energy range of interest labelled. Right'. ACX model folded through the Chandra response, 
with the main charge exchange lines labelled. In this model the relative abundances of all of the metals are set equal to one another. 



Energy (keV) Energy (keV) 

Figure 7. Top left: Same spectrum as in Fig. [4] but with a charge exchange model overplotted with different levels of absorption, allowing us to compare 
the spectral shape. The solid black line is for a model with the LAB survey level of absorption (0.14xl0 22 cm -2 ). Letting the absorbing column be free 
and fitting results in the solid red line, for which the column is 0.35xlO 22 cm -2 . Top right: Fitting the same spectrum but now adding an absorbed thermal 
component (dashed black curve) to the charge exchange model (dashed blue curve). The fit remains very poor. Bottom left: Fitting the filament spectrum with 
an absorbed variable abundance charge exchange model. The colour scheme of the lines is the same as in Fig. [6] Bottom right: Fitting with an absorbed variable 
abundance charge exchange model (blue dashed curve) plus an absorbed apec model (black dashed curve). 


To try to improve the fit further, particularly above 1 keV, we 
add an absorbed APEC component, so that the model becomes 
PHABS(VACX+APEC). Fitting for the broad 0.5-7.0 keV band, 
we find that the fit statistic can be improved to a reduced \ 2 of 
1.26 for 446 degrees of freedom, and the best fit spectrum below 


3 keV is shown in the bottom right hand panel of Fig. [7] In this 
fit the relative abundances of 0:Ne:Mg:Si:Fe in the VACX model 
(blue dashed curve) are 0.2:1.4:0:0.4:2.0, with the oxygen abun¬ 
dance again being suppressed relative to Fe. The column density 
is fixed to the LAB survey value of 0.14x 10 22 cm -2 . The APEC 
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component (black dashed curve) is able to improve the fit dramat¬ 
ically above 1 keV. The temperature of the APEC component is 
2 . 2 lg '2 keV, which is reasonable. In this model, the charge ex¬ 
change component is the dominant component in the 0.5-1.0keV 
band, contributing 57 percent of the X-ray flux. 

If charge exchange is a significant contributor to the 0.5-1.0 
keV band, the constant low temperature, ~0.8 keV, component seen 
in the two thermal component model fits to the filament spectra 
(Fig. CD may be modelling line in the charge exchange spectrum 
originating from the Lyman /?, 7 and <5 transitions in OVIII. This 
would provide an explanation for why the same low temperature 
component is always found in the two temperature thermal fits to 
the spectra. 

We therefore conclude that a model consisting of charge ex¬ 
change plus an APEC model is a possibility for explaining the 
Chandra ACIS filament spectrum of the Perseus cluster, with 
charge exchange being the dominant contributor to the X-ray flux 
in the 0.5-1.0 keV band. While the fit is not formally as good as for 
an absorbed two temperature thermal model, further improvements 
in the charge exchange model may result in an improved fit. Further 
theoretical work is also needed to determine whether the normal¬ 
isation of the charge exchange emission can match the observed 
X-ray flux from the filaments I Fabian et al.l201 lh . 

We note that lWemer et all d2013l) have also previously exam¬ 
ined the soft X-ray emission from the filaments in the Virgo cluster, 
and found that they are brightest in the 0.7-0.9 keV band. This is 
precisely the band in which the ACX charge exchange model pre¬ 
dicts there to be significant contributions from the Ly-/3, 7 and S 
lines of OVIII. 


4.2 Comparing XMM RGS spectra of Centaurus filaments 

Here we compare the RGS spectrum of the X-ray emitting filamen¬ 
tary plume in the Centaurus cluster with the ACX charge exchange 
model. In the XMM observations we study, (obsids:0046340101 
and 0406200101), the roll angle has been set to make the plume 
direction perpendicular to the RGS dispersion direction (shown in 
the left hand panel of Fig. [ 8 j>. This arrangement allows an RGS 
spectrum for the plume to be obtained which excludes contribu¬ 
tions from the nucleus, and this spectrum is shown in the right 
hand panel of Fig. [ 8 ] The data reduction of the RGS data followed 
the methods described in lSanders et al.1120081) . The coordinates of 
the extraction position for the spectrum are RA,DEC= (192.2194,- 
41.303982), and we used an inner 90 per cent PSF extraction, cor¬ 
responding to a 60 arcsec wide strip across the cluster. 

In the right hand panel of Fig. [ 8 ] we compare the RGS spec¬ 
trum of the plume with the ACX charge exchange spectrum. The 
plume spectrum does not appear to contain obvious evidence of 
emission from the OVII line which would be expected to be promi¬ 
nent. The Ly-/?, Ly -7 and Ly-<5 lines from OVIII are also not clearly 
present. This suggests that charge exchange is not the dominant 
contributor to the X-ray emission from the filamentary plume in 
the Centaurus cluster. 

We emphasise, however, that this X-ray plume does not con¬ 
tain any H-a filaments, as can be seen by comparing the X-ray 
image in Fig. [ 8 ] with the H-a image of Centaurus in the bottom left 
hand panel of Fig. [7] The RGS spectrum for Centaurus we present 
here is therefore not extracted from the same type of environment 
as the Chandra spectrum of the Perseus filaments we presented in 
section |4~T1 This may provide an explanation for the weakness of 
the any charge exchange lines in the RGS spectrum. We also stress 
that the RGS spectrum has a significant continuum contribution 


from the surrounding intracluster medium, which cannot be sub¬ 
tracted accurately, and which may obscure any emission lines from 
charge exchange. The effective area of the RGS is also reduced in 
the 20-24 A region compared to the rest of the spectrum due to the 
malfunction of CCD4 of the RGS2 detector, lowering the signal to 
noise. 


5 ABSORPTION BY COLD GAS 

To explore whether the gas in the filamentary structure be seen in 
absorption in X-rays, we explore a region in the Perseus cluster, 
shown by the white box in Fig. [9] This region has been the focus 
of detailed CO observations I lSalome et alJl200S) . which found an 
average gas mass in this region of ~5 x 10 7 Mq. Taking the volume 
of the molecular gas region to be a cuboid of dimensions lOOOpc 
x 500pc x 500 pc, this yields a number density of the molecular 
gas of ~10 particles cnC 3 . Taking the depth to be the same as its 
width (500pc), this yields a sizeable expected column density of 
10 22 cm -2 , which we would expect to be able to discern from X- 
ray observations. 

To test whether this absorption is seen in the X-ray obser¬ 
vations, we compared the X-ray spectrum on the filament region 
(the white box in Fig. [9} with a neighbouring region devoid of 
filaments (the yellow box in Fig. [9}, which also lies outside the 
X-ray cavities in a region with similar X-ray surface brightness. 
These spectra are shown in Fig.|T0] where the red spectrum is from 
the region with molecular gas we study (the white box in Fig. ED. 
and the black spectrum is for the nearby offset region (the yellow 
box in Fig. ED. Both spectra have been scaled to account for differ¬ 
ences in their extraction areas. The best fit model shown is fit to the 
offset region only, and the lower panel of Fig. [TO] shows the ratio 
to this model. This model is an absorbed two component thermal 
model, PHABS(APEC1+APEC2), where the column density in the 
PHABS component is fixed to the LAB survey value of 0.14x 10 22 
cm -2 . We see that there is clear deficit in the X-ray spectrum of the 
molecular gas region below 2 keV, and the spectra are essentially 
identical above 2 keV. 

Qualitatively, this soft X-ray decrement agrees with the idea 
that absorption by the gas in the filament region is taking place. 
However the spectral shape of the decrement cannot be modelled 
accurately using an additional absorption model. As can be seen in 
the ratio plot in Fig. [H the soft X-ray decrement fluctuates, being 
at its greatest in the energy range 0.6-0.7 keV and 1-2 keV, but be¬ 
ing very small in the energy range 0.8-1.0 keV. This variation in the 
decrement below 2 keV cannot be fit for by just adding additional 
absorbing components to the best fit spectrum to the offset region. 

To show that this is the case, we take the best fit model 
to the offset region (the black spectrum in Fig. llOt and attempt 
to fit it to the molecular gas region just by adding a partial 
covering absorption component, so that the model becomes PC- 
FABS *PH AB S (APEC 1 + APEC2). All of the original components 
of the model (the PHABS(APEC1+APEC2) part) are fixed to their 
values for the offset region. The only free parameters are the col¬ 
umn density and covering fraction of the partial covering absorp¬ 
tion model, PCFABS. The result is shown in the top figure of Fig. 
mi (the red spectrum). We see that this model is a poor fit to the 
data. 

One possibility is that in addition to absorption, the fila¬ 
mentary material in the region also emits in X-rays, and this 
acts to partially compensate for the decrement. To test for this 
possibility, we add into the model an additional soft thermal 
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Figure 8. Left : X M M-Newton EPIC image of the Centaurus cluster core, showing that the plume direction is perpendicular to the dispersion direction of the 
RGS. The 90 percent PSF extraction region (a strip of width 60 arcsec) is shown by the dashed white lines. Right'. RGS spectrum of the plume compared to the 
emission lines in the ACX charge exchange spectrum. Note that in this figure the model spectrum has not been broadened to account for spatial broadening, 
and we use it only to show the expected locations of the lines. The spectrum also includes significant continuum emission from the surrounding intracluster 
medium. We see that the prominent OVII line and the OVIII I .y-y and Ly-<5 lines expected from charge exchange are not present. 


component to represent possible filament emission. The model 
is now PCFABS *PHABS( APEC 1+APEC2+APEC3), where the 
PHABS(APEC1+APEC2) is fixed to the best fitting values to the 
offset spectrum as before. The free parameters are the column den¬ 
sity and covering fraction of the partial covering PCFABS model, 
and the temperature and normalisation of the new thermal compo¬ 
nent (APEC3) describing possible filament emission. 

This model produces a good fit to the data (a reduced \ 2 
of 1.07 for 684 degrees of freedom), and is shown in the bot¬ 
tom panel of Fig. El The best fitting component values are very 
realistic: the partial covering model gives a column density of 
2 ±° | x 10 22 cm” 2 (similar to the ~ 10 22 cm” 2 we estimated 
earlier from the mass of molecular gas present) with a covering 
fraction of 0.15lg o 2 . The temperature of the added APEC com¬ 
ponent is l.Olg'gg keV (with a metal abundance of 0.5 Zq) which 
is reasonably similar to what is obtained by fitting a typical fila¬ 
ment spectrum from Perseus with a single absorbed APEC com¬ 
ponent. The intensity of this component is 4.3 xlO” 16 ergs cm” 2 
s” 1 arcsec” 2 in the 0.6-2.0 keV band, which is very similar to the 
average intensity of the 12 filaments in Perseus we studied earlier, 
which is 4.0x 10” 16 ergs cm” 2 s” 1 arcsec” 2 in the same band. 

We therefore conclude that the X-ray spectrum of the molecu¬ 
lar region is consistent with absorption from molecular gas, which 
is partially compensated for by X-ray emission from the filaments. 


6 CONCLUSIONS 

We have explored the X-ray properties of the filaments surrounding 
the BCGs in four nearby, X-ray bright galaxy clusters using deep 
Chandra observations: the Perseus cluster, the Centaurus cluster, 
the Virgo cluster and Abell 1795. In the Perseus cluster, we have 
studied the X-ray properties of 12 individual filaments, all of which 
can be fitted with an absorbed two component thermal model. The 



Figure 9. Composite image of the centre of the Perseus cluster. The Ha 
filaments are shown in pink, while the X-ray emission from Chandra obser¬ 
vations is shown in green. The white box is located on a prominent area of 
molecular gas as studied by lSalome et alj <2008 ). and we compare the X- 
ray spectrum from this region with that from a nearby region (yellow box) 
with no intervening molecular gas along the line of sight. 


temperatures of these two components are consistent within errors 
between all of the 12 filaments, with a lower temperature of 0.77 
keV, and an upper temperature 1.77 keV. This absorbed two com- 
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Figure 10. Comparing the X-ray spectra of the molecular gas region (red 
spectrum, white box in Fig. [9} with that of an offset region (black spectrum, 
yellow box in Fig. [9j. Both spectra have been scaled to account for any 
differences in the areas in the extraction regions. The best fit model is fit 
to the black spectrum (for the offset region), and the lower panel shows the 
residuals to this model. It is clear that in the region with molecular gas (red 
spectrum) there is a deficit in the X-ray emission at energies below 2 keV. 



Figure 11. The same spectra as shown in Fig. 1101 but focusing just at the 
0.6-2.0 keV soft band. Top: Here we attempt to fit the spectrum for the 
molecular region (red) by just adding a partial covering absorption model 
to the best fit spectrum to the offset region (the black spectrum). As the 
residuals show, the resulting fit is very poor. Bottom: Same as the top panel 
but we now add an additional APEC component to represent possible X- 
ray emission from the filaments (shown as the lowest dashed red line in the 
spectral model). This results in a good fit. 


ponent thermal model also proved to be a good fit to the spectra of 
the stacked filaments in Virgo, Centaurus and Abell 1795, and sim¬ 
ilar temperatures of the two components were found. The filaments 
in all of these systems therefore appear to have the same spectral 
shape, suggesting that the physical process causing the soft X-ray 
emission from the filaments is the same in all of these systems. 

We find that it is also possible to fit the stacked Chandra ACIS 
filament spectrum from Perseus with a model consisting of a charge 
exchange component and a soft thermal component. To obtain a 
good fit below 1 keV, the oxygen abundance needs to be suppressed 
relative to Fe, to reduce the intensity of the OVII and OVIII lines, an 
effect which cannot be achieved by just raising the level of absorp¬ 
tion. Below 1 keV, it is possible to adequately fit the filament spec¬ 
trum with just a variable abundance charge exchange model (bot¬ 
tom left hand panel of Fig.[7J. Flowever, to obtain the correct spec¬ 
tral shape above IkeV, an additional low temperature APEC com¬ 
ponent is needed (bottom right hand panel of Fig. [7}. The charge 
exchange still provides the dominant contribution to the spectrum 
in the 0.5-1.0 keV band. Whilst the fit quality is not as good as an 
absorbed two temperature thermal model, further improvements to 
the charge exchange model in the future may improve the quality of 
the fit. These findings therefore add support to the idea that charge 
exchange at least contributes to the X-ray emission we see from the 
filaments. 

We find that when we examine the high spectral resolution 
RGS spectrum of the filamentary plume seen in X-rays in the Cen¬ 
taurus cluster, the OVII line, which is expected to be prominent in 
the charge exchange spectrum, is not present. It therefore appears 
unlikely that the soft X-ray emission seen in the X-ray plume in 
Centaurus is dominated by charge exchange emission. This X-ray 
plume does not however contain any H-o? filaments, making it fun¬ 
damentally different to the filamentary regions of the Perseus clus¬ 
ter we explore, and possibly explaining the weakness of any charge 
exchange spectral features. The RGS spectrum also contains a sig¬ 
nificant background contribution from the surrounding ICM which 
cannot be accurately removed, and may prevent us from detecting 
the OVII line. 

In the Perseus cluster, we search for signs of X-ray absorption 
due to gas in the filaments, using a region which has a large amount 
of molecular gas, and comparing it to a nearby offset region where 
there are no filaments. At energies below 2 keV, there is a deficit 
in the X-ray spectrum of the molecular gas region, suggesting that 
absorption is playing a role. However the X-ray spectrum of this 
deficit is complex, and cannot be fit by just simply adding addi¬ 
tional absorbing components. 

This suggests that, in addition to absorption, the gas also emits 
X-rays which partially compensates for the absorption, leading to 
a more complex spectral shape. We find that when we add an addi¬ 
tional soft thermal component to represent X-ray filament emission, 
we obtain a good fit to the X-ray spectrum of the molecular gas re¬ 
gion. Moreover, the physical parameters of the best fit model are 
very realistic. The column density of the absorption model agrees 
well with what would b e expected giv en the molecular gas mas 
in the region obtained by lSalome et alj i2008l) . The X-ray surface 
brightness of the filament component is in good agreement with 
the average X-ray surface brightness of the 12 filaments we have 
studied in Perseus. 
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Table Al. Chandra data used in this paper. 


Object 

Obs ID 

Exposure (ks) 

RA 

Dec 

Start Date 

Perseus 

3209 

95.77 

03 19 47.60 

+41 30 37.00 

2002-08-08 


4289 

95.41 

03 19 47.60 

+41 30 37.00 

2002-08-10 


4946 

23.66 

03 19 48.20 

+41 30 42.20 

2004-10-06 


4947 

29.79 

03 19 48.20 

+41 30 42.20 

2004-10-11 


6139 

56.43 

03 19 48.20 

+41 30 42.20 

2004-10-04 


6145 

85.00 

03 19 48.20 

+41 30 42.20 

2004-10-19 


4948 

118.61 

03 19 48.20 

+41 30 42.20 

2004-10-09 


4949 

29.38 

03 19 48.20 

+41 30 42.20 

2004-10-12 


6146 

47.13 

03 19 48.20 

+41 30 42.20 

2004-10-20 


4950 

96.92 

03 19 48.20 

+41 30 42.20 

2004-10-12 


4951 

96.12 

03 19 48.20 

+41 30 42.20 

2004-10-17 


4952 

164.24 

03 19 48.20 

+41 30 42.20 

2004-10-14 


4953 

30.08 

03 19 48.20 

+41 30 42.20 

2004-10-18 

Centaurus 

504 

31.75 

12 48 48.70 

-41 18 44.00 

2000-05-22 


4954 

89.05 

12 48 48.90 

-41 18 44.40 

2004-04-01 


4955 

44.68 

12 48 48.90 

-41 18 44.40 

2004-04-02 


5310 

49.33 

12 48 48.90 

-41 18 44.40 

2004-04-04 


16223 

180.0 

12 48 48.90 

-41 18 43.80 

2014-05-26 


16224 

42.29 

12 48 48.90 

-41 18 43.80 

2014-04-09 


16225 

30.1 

12 48 48.90 

-41 18 43.80 

2014-04-26 


16534 

55.44 

12 48 48.90 

-41 18 43.80 

2014-06-05 


16607 

45.67 

12 48 48.90 

-41 18 43.80 

2014-04-12 


16608 

34.11 

12 48 48.90 

-41 18 43.80 

2014-04-07 


16609 

82.33 

12 48 48.90 

-41 18 43.80 

2014-05-04 


16610 

17.34 

12 48 48.90 

-41 18 43.80 

2014-04-27 

Virgo 

5826 

126.76 

12 30 49.50 

+12 23 28.00 

2005-03-03 


5827 

156.20 

12 30 49.50 

+12 23 28.00 

2005-05-05 


5828 

32.99 

12 30 49.50 

+12 23 28.00 

2005-11-17 


6186 

51.55 

12 30 49.50 

+12 23 28.00 

2005-01-31 


7210 

30.71 

12 30 49.50 

+12 23 28.00 

2005-11-16 


7211 

16.62 

12 30 49.50 

+12 23 28.00 

2005-11-16 


7212 

65.25 

12 30 49.50 

+12 23 28.00 

2005-11-14 

Abell 1795 

See tables Al and A2 from Walker et al. (2014) 





Coma 

555 

8.66 

12 59 48.00 

+27 58 00.00 

1999-11-04 


1112 

9.65 

12 59 48.00 

+27 58 00.00 

1999-11-04 


1113 

9.65 

12 59 48.00 

+27 58 00.00 

1999-11-04 


1114 

9.05 

12 59 48.00 

+27 58 00.00 

1999-11-04 


9714 

29.65 

12 59 48.00 

+27 58 00.00 

2008-03-20 


13993 

39.56 

12 59 51.92 

+27 54 17.76 

2012-03-21 


13994 

81.99 

12 59 39.42 

+27 54 28.59 

2012-03-19 


13995 

62.98 

12 59 59.57 

+27 55 02.47 

2012-03-14 


13996 

123.06 

12 59 51.01 

+27 56 33.68 

2012-03-27 


14410 

78.53 

12 59 51.92 

+27 54 17.76 

2012-03-22 


14411 

33.64 

12 59 39.42 

+27 54 28.59 

2012-03-20 


14415 

34.53 

12 59 59.57 

+27 55 02.47 

2012-04-13 


14406 

24.76 

12 59 59.57 

+27 55 02.47 

2012-03-15 

Abell 3627 

8178 

57.41 

16 15 03.80 

-60 54 25.20 

2007-07-08 


9518 

140.04 

16 13 25.59 

-60 45 43.10 

2008-06-13 


12950 

89.85 

16 13 39.28 

-60 52 10.20 

2011-01-10 
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